Chandra Observation of the Cluster of Galaxies MS 0839.9+2938 
at z = 0.194: the Central Excess Iron and SN la Enrichment 



Yu Wang 1 ' 2 ' 3 , Haiguang Xu 1 , Zhongli Zhang 1 , Yueheng Xu 1 , Xiang-Ping Wu 2 , Sui-Jian 

Xue 2 , and Zongwei Li 3 

ABSTRACT 

We present the Chandra study of the intermediately distant (z = 0.194) 
cluster of galaxies MS 0839.9+2938. By performing both the projected and 
deprojected spectral analyses, we find that the gas temperature is approximately 
constant at about 4 keV in 130 — 444/iyQ 1 kpc. In the inner regions, the gas 
temperature descends towards the center, reaching <3 keV in the central 37hj^ 
kpc. This infers that the lower and upper limits of the mass deposit rate are 
9 — 34M yr -1 and 96 — 126M Q yr _1 , respectively within 74/iyQ 1 kpc where the gas 
is significantly colder. Along with the temperature drop, we detect a significant 
inward iron abundance increase from about 0.4 solar in the outer regions to ~1 
solar within the central 37h^Q kpc. Thus MS 0839.9+2938 is the cluster showing 
the most significant central iron excess at z>0.2. We argue that most of the 
excess iron should have been contributed by SNe la. By utilizing the observed 
SN la rate and stellar mass loss rate, we estimate that the time needed to enrich 
the central region with excess iron is 6.4 — 7.9 Gyr, which is similar to those 
found for the nearby clusters. Coinciding with the optical extension of the cD 
galaxy (up to about 30h^ kpc), the observed X-ray surface brightness profile 
exhibits an excess beyond the distribution expected by either the (3 model or the 
NFW model, and can be well fitted with an empirical two-/? model that leads to 
a relatively flatter mass profile in the innermost region. 

Subject headings: galaxies: clusters: individual (MS 0839.9+2938) — intergalac- 
tic medium — galaxies: evolution — X-rays: galaxies: clusters 
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1. Introduction 

The metallicity of the intra-cluster medium (ICM) in clusters of galaxies is a sensitive 
tracer of both the history of the hierarchical clustering and the chemical evolutions of mem- 
ber galaxies and clusters themselves. Since the ICM has been enriched by the supernova 
explosions, the study of the metal abundance gradient in clusters, especially in those distant 
ones, may provide us with crucial constraints on the SN rate, stellar mass loss rate, and thus 
star formation rate. Along with the metal enrichment, a tremendous amount of the super- 
nova energy has been inputted into the ICM, primarily in the form of gas dynamic energy 
that may be used to heat the ICM. Thus detailed measurements of the metal distribution 
in clusters as a function of redshift may also provide us with valuable hints to investigate 
how the ICM has been heated, a puzzling question that has plagued astronomers for a long 
time since the absence of massive cooling gas has been confirmed in many clusters (see, e.g., 
Makishima et al. 2001 for a review). 

In the past five years, the inward iron abundance increase, first discovered with Ginga, 
ROSAT and ASCA in some nearby clusters that usually host a giant cD galaxy (e.g., Koyama 
et al. 1991; Fukazawa et al. 1994; Ezawa et al. 1997; Xu et al. 1997; Ikebe et al. 1999; 
Finoguenov et al. 2000), has been revealed in more sources (e.g., a sample of 12 nearby rich 
clusters at z<0.1 that have cool cores; De Grandi et al. 2004) and in much more details by 
the high spatial resolution observations of Chandra and XMM, which allows us to measure 
the abundances of iron and other metals with sufficient accuracy to deduce the SN la and 
SN II enrichments and their spatial variations. In clusters that exhibit a central iron excess, 
it has been found that the average iron abundance is typically 0.4 — 0.6 solar within a 
radius of 30-80 kpc, which descends to about 0.2 — 0.3 solar in the outer regions. In several 
most neighboring sources, such as M87/the Virgo cluster (Matsushita et al. 2003), NGC 
1399/the Fornax cluster (Buote 2002), NGC 5044/the NGC 5044 group (Buote et al. 2003), 
NGC 1129/AWM 7 (Furusho et al. 2003), and NGC 4696/the Centaurus cluster (Sanders 
& Fabian 2002), the average iron abundance is found to increase further to 0.7-1.5 solar 
within the central 15-35 kpc, although in the innermost 2 kpc of M87 and NGC 4696 the 
iron abundance decreases down to <0.5 solar to form a central dip. In M87 (Matsushita et al. 
2003), NGC 1399 (Buote 2002), NGC 5044 (Buote et al. 2003) and a sample of 19 clusters 
(Tamura et al. 2004), the spatial distributions of the silicon and sulphur abundances are 
likely to follow that of iron. The distribution of the oxygen abundance, on the other hand, 
is roughly uniform within the errors. These results strongly support the idea that there has 
been an excess SN la enrichment in the cluster's inner region. 

For the time being, the spatial distributions of metal abundances in clusters at z > 0.1 
are less understood due to the instrumental limitations, although there is evidence that the 
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same central iron abundance excess may also exist (e.g., Allen et al. 2001, 2002; Iwasawa et 
al. 2001). As a step to approach more distant clusters, in this paper we present the Chandra 
study of the iron abundance distribution in MS 0839.9+2938, an intermediately distant 
cluster located at z — 0.194 (Stocke et al. 1991). The cluster was listed as an extended 
X-ray source and identified as a cluster in the Einstein Extended Medium Sensitivity Survey 
(Gioia et al. 1990). It hosts a giant cD galaxy at its dynamic center, which is a weak, 
compact radio source with an intense, extended H a emission (Nesci et al. 1995). The advent 
of studying this source is that since its average gas temperature is low (~ 3.6 keV; see §4.2), 
it is easier to measure the abundance gradient accurately due to the large equivalent widths 
of Fe-L and -K lines. 

We organize the paper as follows. In §2, we describe the observation and data reduction. 
In §3 and §4, we present the imaging and spectral analyses, respectively. In §5, we investigate 
the spatial distributions of the gas and gravitating mass as well as gas cooling. Finally, we 
discuss and summarize our results in §6 and §7, respectively. Throughout the paper, we 
adopt the cosmological parameters Ho = 70 km s _1 Mpc -1 , Q m = 0.3 and Q\ = 0.7, so that 
1" corresponds to ~ 3.2/iyQ 1 kpc at the redshift of the cluster. In order to compare directly 
with previous results, we use the older solar abundances from Anders & Grevesse (1989), 
where the iron abundance relative to hydrogen is 4.68 x 10 -5 in number, which is about 46% 
higher than the currently preferred meteoritic value (McWilliam 1997; Grevesse & Sauval 
1999). Unless stated otherwise, the quoted errors are the 90% confidence limits. 



2. Observation and Data Reduction 

The Chandra observation of MS 0839.9+2938 was carried out on January 29, 2001 for 
a total exposure of 29.7 ks with the chips 2, 3, 6, 7, 8, and 9 of the Advanced CCD Imaging 
Spectrometer (ACIS). The center of the cluster was positioned close to the nominal aim 
point on the S3 chip (CCD 7), so that most of the X-ray emission of the galaxy was covered 
by the S3 chip. The events were collected with frame times of 3.2 s and telemetered in the 
Faint mode. The focal plane temperature was set to be -120 °C. In this work, we used the 
Chandra data analysis package CIAO software (version 3.0) to process the data extracted 
from the S3 chip only. We kept events with ASCA grades 0, 2, 3, 4, and 6, and removed all 
the bad pixels, bad columns, and columns adjacent to bad columns and node boundaries. 
In order to detect occasional background flares, whose effects are particularly significant on 
the backside-illuminated SI and S3 chips, we examined the lightcurves of the background 
regions on the S3 chip both in 2.5-7 keV, where the background flares are expected most 
significant, and in 4-7 keV, where the source contamination is least. We found that there was 
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almost no strong background flares that increased the count rate to over 20% more than its 
mean value. By excluding the high background intervals we obtained a net exposure of 27.1 
ks for the analyses. We restrict our spectral study in 0.7-8.0 keV to avoid the calibration 
uncertainties at the low energy end and the instrumental background at the high energy 
end. In the spectral fittings, we have taken into account corrections for the ACIS quantum 
efficiency degradation and used the XSPEC vll.2.0 software. 



3. Imaging Analysis 

3.1. X-ray Morphology 

In Figure 1, we present both the smoothed 0.7-8.0 keV image of the central 4.9' (944/ifQ 1 
kpc) of MS 0839.9+2938 in the logarithmic scale, and the corresponding DSS optical image on 
which the X-ray contours are overlaid. The smoothed X-ray image was created by using the 
CIAO tool csmooth with a minimum significance of 3 and a maximum significance of 5, and 
has been background-subtracted and exposure-corrected. As can be seen clearly, the diffuse 
X-ray emission is strongly peaked at the cluster's center. Within the innermost 8" (~ 26/i^q 1 
kpc), it appears to be elongated approximately in the northwest-southeast direction, where 
it follows the distribution of the optical lights. Outside this region the diffuse emission is 
elongated in the north-south direction out to 1.1' (~ 210h^ kpc). The X-ray emission of the 
cluster extends to at least 4.2' (~ 800/ifQ 1 kpc), where the background begins to dominate 
the flux. There is no apparent substructures in X-rays, which indicates that the cluster has 
not experienced any major merger events in the recent past, and is likely to be in a relaxed 
state. Since the astrometric errors of the Chandra observation are about 1", the position 
of the X-ray peak (RA=08h42m56.0s, Dec=+29d27m27.2s, J2000) is in excellent agreement 
with the optical center of the cD galaxy 2MASX J08425596+2927272 (Becker et al. 1995) 
to within < 0.5". 

Based on the studies of Brandt et al. (2000) and Mushotzky et al. (2000), we estimate 
that there may be 5-10 unrelated X-ray sources in the r < 4.2' region of the cluster. By 
utilizing a wavedetect detection algorithm and a source detection threshold of 10~ 6 , and then 
crosschecking the results both by using the celldetect detection algorithm and by eye, we 
detected 10 X-ray point sources in r < 4.2', among which less than 1 is expected to be a false 
detection. We find that two of them have optical and radio counterparts within 1". One is 
CGCG 150-019, an Sb galaxy located at z = 0.0276 (Wegner et al. 2001), and the other 
is FIRST J084251. 9+292825, a Seyfert galaxy in MS 0839.9+2938 (Hutchings & Edwards 
2000). In the analysis that follows, all the detected X-ray point sources were excluded. 
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3.2. Central Excess Emission 

In Figure 2, we show the azimuthally-averaged surface brightness profile of the cluster 
in 0.7-8 keV that has been corrected with a weighted exposure map. All the detected point 
sources have been excluded, while the background has not been subtracted. We find that 
in 15.8 — 800/iyQ 1 kpc, the observed X-ray surface brightness profile can be well described 
(xt — 1-05) by a model consisting of a standard (3 component (Jones & Forman 1984) and 
a spatially uniform background as 

S(r) = S [l + (r/r c ) 2 } ^ + S hkg , (1) 

where r c is the core radius, (5 is the slope, Sbk g is the background, and So is the normalization. 
The best-fitting parameters are r c = 51.7 ± 0.5 kpc and f3 = 0.62 ± 0.01. However, 
when we extrapolate the model inwards, the fit becomes worse (r c = 41.7 ± 0.4 /ifo 1 kpc, 
(3 = 0.59 ± 0.01, and xt = 2-84), which significantly underestimates the data in r < 4.9" 
(15.8/ifo 1 kpc), leaving an obvious central excess emission above the model. Such central 
excess emission has been seen in many clusters, e.g., the Centaurus cluster (Ikebe et al. 
1999), Abell 1795 (Xu et al. 1998), Abell 1983 (Pratt & Arnaud 2003), and C10024+17 
(Ota et al. 2004). It may be attributed to either a temperature drop, or a metal abundance 
increase, or an excess of the gas density, all within the central region. In §5.1, we will present 
detailed model fittings to describe the observed X-ray surface brightness profile by taking into 
account both the temperature and abundance gradients obtained from the spectral analyses 
(§4). 

4. Spectral Analysis 

4.1. Two-Dimensional Hardness Ratio Distribution 

Before the direct spectral fittings, we first study the two-dimensional hardness ratio 
distribution of the hot diffuse gas in the cluster. The hardness ratio HR is defined as the 
ratio of the background-subtracted and exposure-corrected counts in 1.4-8.0 keV to those 
in 0.7-1.4 keV. As can be seen in Figure 3, in ~ 100 — 300/ifo 1 kpc the hardness ratio is 
approximately constant at a value of 1.00 — 1.30, with typical errors of ±0.06. Assuming an 
average metallicity of 0.4 solar and the Galactic absorption A" H = 4.04 x 10 20 cm~ 2 (Dickey 
& Lockman 1990), this is consistent with the hardness ratio of an isothermal plasma at kT ~ 
4.0 — 6.5 keV. In r ~ 70 — 140/ifo kpc, there are two relatively colder (90% confidence level) 
regions located at the southeast and northwest of the center, respectively. In the southeast 
colder region the gas temperature is estimated to be 2.8 — 3.2 keV (HR= 0.80 ± 0.04), which 
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is close to that of the gas in the northwest colder region (3.0 — 3.2 keV, HR= 0.83 ± 0.04). 
In r<40hjQ kpc, where the hardness ratio is ~ 0.60 — 0.83 with typical errors of ±0.04, the 
gas appears to be significantly colder than in the outer regions, indicating the existence of 
a plasma at kT ~ 2.2 — 3.1 keV. These results agree very well with those obtained in the 
detailed spectral fittings (§4.3). 

4.2. Total Spectrum 

We extracted the total spectrum of the cluster using a circular region filter with a radius 
of ~ 2.3' (~ 440/itq 1 kpc), which is centered at the X-ray peak, excluding all the detected 
point sources. The background spectrum was extracted from a separate region on the S3 
chip that is far away enough from the cluster. We also have used the Chandra blank fields 
for background, and obtained essentially the same results as those shown below. At the 
boundary of the so defined cluster region, the count flux is about 3 times more than that of 
the background. The resulting background-subtracted spectrum contains a total of about 
14,100 counts in 0.7-8.0 keV, showing a strong line feature at ~ 5.6 keV that corresponds to 
a rest-frame energy of about 6.7 keV where the K a lines of He-like iron should reside. The 
appearance of this Fe K a feature indicates that the average gas temperature of the cluster 
should be within 2.7-5.2 keV. 

We have fit the spectrum with a model consisting of a single absorbed isothermal com- 
ponent (Table 1). First we choose the APEC model (Smith et al. 2001) as the isothermal 
component, and fix the redshift and absorption to z = 0.194 and the Galactic value, respec- 
tively. The fit is marginally acceptable (xt = 1-18), which gives an average gas temperature 
of kT = 3.65^Q2o keV, and an average metal abundance of Z = 0.44^QQg solar. In order to 
improve the fit we have attempted to apply the VAPEC model instead of the APEC model. 
For MS 0839.9+2938, the redshift is large enough to move the Ly a lines of silicon to lower 
energies to avoid the calibration uncertainties caused by the mirror Ir edge. Thus in the 
fitting, we leave both the Fe and Si abundances free, and tie the abundance of Ni to that 
of Fe. Although we have applied the contamination file from the latest CALDB to correct 
for the ACIS deficiency in low energies, the degradation of energy resolution still does not 
recover. For an approximate estimation, we tie the abundances of other a-elements (O, Ne, 
Mg, S, Al, Ar and Ca ) together. The model yields a gas temperature of kT = 3.56jlo'|o keV, 
and abundances of Zp c = 0.46^QQg solar, Z$i = 0.42+Q2g solar, and Z a = 1.201q42 s °l ar that 
is mostly dominated by the S abundance. Although the VAPEC fit (xl — 1-08) is better 
than the APEC fit, there are still residuals remaining at about 0.9 keV, where the model 
underestimates the data. These residuals may mostly caused by the excess Fe-L emissions 
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arising from the colder gas in the central region, whose existence has been inferred in the 
study of the hardness ratio distribution (§4.1). In either the APEC or VAPEC fitting, even 
if the absorption is left free, the obtained column density is still consistent with the Galactic 
value, and the fit is not improved. 

We also have examined if the fit to the total spectrum can be further improved by 
the use of a model consisting of two isothermal spectral components, both subjected to a 
common absorption (Table 1). For this purpose we employ two VAPEC components, whose 
absorption is again fixed to the Galactic value because allowing it to vary does not improve 
the fit. The metal abundances of each VAPEC component are grouped in the same way as 
in the single VAPEC fitting, and the abundances of the same element of the two VAPEC 
components are tied together. We perform the F-test and find that at the 95% confidence 
level the two-temperature fit (xl — 1-01) is better than the one-temperature APEC or 
VAPEC fit. The obtained gas temperatures for the two components are kT\ ow = 2.10^°;g4 
keV and A;T high = 6.041^24 keV, and the metal abundances are Z Fc = 0.38^^1 solar, Z Si = 
0.371q24 solar, and Z a = 0.96lo' 5O solar. With these best-fit parameters we calculate that 
the unabsorbed 0.5-10 keV flux of the cluster is 3.9 jloll x 10~ 12 erg cm -2 s^ 1 , corresponding 
to a luminosity of 4.1+^2 x 10 44 /if 2 erg s _1 . 



4.3. Projected and Deprojected Analyses 

In order to investigate the temperature and abundance gradients, we divided the cluster 
into five annular regions, all centered at the X-ray peak. As the first step, we study the 
projected spectra extracted in these annuli. The background spectrum was again extracted 
in a region on the S3 chip far away enough from the cluster; the use of the Chandra blank fields 
yields essentially the same results. Assuming that the gas in each annulus is in the thermal 
equilibrium, we fit the spectra using a single absorbed APEC model with the absorption 
fixed to the Galactic value; allowing the absorption to vary does not improve the fits. As is 
shown in Figure 4 and Table 2, the fit is acceptable for each region. The gas temperature 
is found to be approximately constant at ~ 4.3 keV outside about 130/ifg 1 kpc, while it 
descends to ~ 3.2 keV in the central 37/ifg 1 kpc, as has been expected in the study of the 
hardness ratio distribution (§4.1). On the other hand, the average metal abundance, which 
is mainly determined by iron, increases towards the center from a roughly constant value of 
0.4 solar in the outer regions to about 0.8 solar within the central 37/iyQ 1 kpc. Instead of 
the APEC model, we have also attempted to apply the VAPEC model (Fig. 4 and Table 
2). For better statistics, all the abundances of the a-elements are tied together, and the 
absorption is again fixed to the Galactic value. We find that the fits are acceptable too, and 
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the obtained gas temperatures agree nicely with those obtained with the APEC model. In 
the regions outside r = 74/iyg 1 kpc, the iron abundance varies slowly in a range of 0.25-0.42 
solar, while in the central region, it increases drastically to about 0.8 solar. The spatial 
distribution of the average abundance of the a-elements, however, is not well determined 
due to large errors. 

Next, we perform the deprojected analysis of the spectra extracted in the inner four 
regions by using the standard "onion peeling" method (Fig. 4 and Table 2). Assuming 
that the gas temperature and abundance are uniform in each spherical shell, and that the 
gas density follows the best-fit profile derived in §5.1, we subtract emissions from the outer 
shells projected on the inner shells. We fit the deprojected spectra with either the APEC or 
VAPEC model whose absorption is fixed to the Galactic value. We find that the deprojected 
spectra can be well fitted with either of the models. The gas temperatures obtained with 
the APEC model agree very well with those obtained with the VAPEC model, and both the 
average metal abundance of the APEC model and the iron abundance of the VAPEC model 
show significant central increase in the innermost region as has been found in the projected 
analysis. For each model the best-fit gas temperatures and metal abundances obtained in 
the projected and deprojected analyses are consistent with each other. In the innermost 
region, the deprojected temperatures are only slightly lower than, but still consistent with 
the projected values. 



4.4. Temperature and Iron Abundance Gradients in the ICM 

As we have shown above, both the projected and deprojected spectral analyses indicate 
that in MS 0839.9+2938 there is a central temperature drop within a radius of about 74h^Q 
kpc. According to the deprojected APEC fit, the gas temperature descends from kT ~ 4.2 
keV in > 130/iyQ 1 kpc to ~ 2.7 keV in < 37hjQ kpc, which can be approximated by a power- 
law profile T{R) = T (R/Ro) a , where T = 3.0 keV, a = 0.17 and R = 35 Ah^ kpc (Fig. 
4a). In the same region where the cold gas is present, the iron abundance shows a clear 
tendency to increase. In the projected analysis, both the iron-dominated metal abundance 
of the APEC model and the iron abundance of the VAPEC model in the central 37hjQ 
kpc region are significantly higher than their counterparts in 74 — 130/ifQ 1 kpc at the 90% 
confidence level. They are also apparently higher than the abundances in 37 — 74/iyQ 1 kpc 
and 130 — 230/iyQ 1 kpc, respectively, with the 90% confidence ranges overlap with each other 
only slightly. In the deprojected analysis, although the obtained errors are much larger, the 
average metal abundance and iron abundance still show a tendency to increase inward. At 
the 90% confidence level, their values inferred for the < 37/ifg 1 kpc region are significantly 
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higher than those for the 74 — l30h 7 Q kpc region. 

In order to crosscheck these results with better statistics, we divide the cluster into the 
inner (< 74/iyQ 1 kpc) and outer (74 — AAAhj^ kpc) regions, and fit their deprojected spectra 
with the an absorbed VAPEC model (Table 3). In the fitting the absorption is fixed to the 
Galactic value, and the metal abundances are grouped in the same way as in the fittings 
of the total spectrum. We find that the fits are acceptable for both regions. At the 90% 
confidence level, the central temperature decrease and abundance increase can be confirmed. 
The silicon abundance and the average abundance of other a-elements, however, show no 
tendency to vary between the two regions. 

Since the iron abundance is determined as much by the Fe-L emission lines keV as by 
the Fe-K lines, the obtained inward iron abundance increase might be an artifact, which is 
actually caused by the enhancement of the Fe-L blend due to an inward temperature drop 
that is even steeper than has been modeled. We have examined this possibility. First we 
exclude the Fe-L lines completely from the spectra and use the Fe-K lines only to measure 
the iron abundance. The obtained iron abundance does show a tendency to increase inwards, 
although it is less significant due to poorer statistics. Second, we fit the deprojected spectrum 
of the innermost region by using a single absorbed VAPEC model with the absorption fixed 
to the Galactic value and the iron abundance fixed to 0.47 solar (the upper limit of the 
iron abundance in the third annular region; Table 2). We do obtain a lower temperature of 
2.38±0.30 keV as is expected. However, the fit appears to be worse and should be ruled out at 
the 90% confidence level in terms of F-test, showing obvious spectral residuals in 4.5 — 7.0 keV 
as well as at ~ 0.9 keV, where the fits to the Fe-K and Fe-L lines are worse. This indicates 
that a central abundance increase is much more preferred than an even steeper temperature 
gradient. Third, we study the two-dimensional fit-statistic contours of temperature and iron 
abundance at the 68%, 90% and 99% confidence levels for the innermost region and the third 
region, respectively, all obtained in the deprojected VAPEC fitting (Fig. 5). As can be clearly 
seen the central iron abundance is higher than that of the third region at a significance of 
99%. Consequently, we conclude that the possibility of the detected central iron abundance 
increase being an artifact should be excluded. Note that the coexistence of the temperature 
drop and abundance increase may indicate that they are physically correlated, since the 
higher the metallicity is, the faster the gas cools down. 
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5. Mass Distributions and Cooling of Gas 

5.1. Gas Density Profile 

Assuming that 1) the gas is ideal and in the thermal equilibrium state, and 2) the 
spatial distribution of the gas is spherically asymmetric, we deduce the three-dimensional 
gas density distribution n g (R) by fitting the observed radial X-ray surface brightness profile 
S(r), which was extracted in 0.7-8.0 keV, azimuthally-averaged, and exposure-corrected. 
The calculations was performed in a projected way by using 

/OO T3J D 

A(T,Z)n 2 g (R)-^=== + S hkg , (2) 

where A(T, Z) is the cooling function calculated by taking into account both the temperature 
and abundance gradients as is modeled by the best-fit deprojected APEC parameters, and 
Sbkg is the background. 

First we assume that the gas density follows the f3 distribution 

n g {R) = n g , [1 + {R/R c fY W/2 , (3) 

where R c is the core radius and (5 is the slope. We find R c = 47.3 ± 0.3hjQ kpc and 
(5 = 0.61 ± 0.01. The fit, however, is unacceptable (xl — 1-63), showing strong residuals in 
< 15.8/ifo kpc where the model significantly underestimates the data (Fig. 2). In fact, the 
gas density derived in such a way is lower than that obtained in the deprojected spectral 
analysis (§4.3) by about 45%. This confirms the result in §3.2, although the central brightness 
excess above the model is less prominent since the effects of the central temperature drop and 
abundance increase are both counted in, which enhance the model emission in the central 
region. 

Next we examine if S(r) can be reproduced by adopting the universal gravitating mass 
distribution as has been found in the cosmological N-body simulations (Navarro et al. (1997) 

p(R) = {R/R s )ll + R/R s r (4) 

where p(R) is the gravitating mass density, R s is the scale radius, p cr u is the critical density 
of the universe at the observed redshift, and 5 C is the characteristic density contract that can 
be expressed in term of the equivalent concentration parameter c as 

200 c 3 
6c = T[ln(l + c)-c/(l + c)r (5) 
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As is shown in Makino et al. (1998) and Wu et al. (2000), under the assumptions of 1) 
spherical symmetry, 2) hydrostatic equilibrium, and 3) ideal and isothermal gas, the gas 
distribution can be given in an analytic form 

n g (x)=n g (0) { V-l ' (6) 

where x = R/R s is the dimensionless radius, a = AnGnm p 8 c p cr i t R 2 s /kT is the index, and 
fi = 0.609 is the average molecular weight for a fully ionized gas. Note that in Eq.(6) a 
background density at infinity has been subtracted to avoid the divergence in integrating 
the X-ray emission along the line of sight. As are shown in Figure 2 and Table 4. We 
find that, comparing with the f3 model, the fit is improved in the central region, giving 
R s = 195.1 ± O.bh^Q kpc, c = 6.5 ± 0.1, and xl — 1-12. However, it still systematically 
underestimates the data in < 15.8 kpc. 

In order to improve the fit in the central region, we have attempted to add an additional 
f3 component to the f3 model to account for the central excess emission, for which the gas 
density distribution is 

n 9 {R) = {n 2 glfi [l + {R/R cl fY m +n% fi [l + {R/ R c2 f] . (7) 

This empirical two-component model has been successfully applied to many clusters of galax- 
ies as well as some groups of galaxies (see, e.g., Makishima et al. 2001 and references therein). 
We find that the model can best fit the observed surface brightness profile throughout the 
spatial range spanned by the cluster and yield an acceptable fit. The best-fit parameters are 
R cl = 65.9+^/ifo 1 kpc and ft = 0.66l^ for one (3 component, which mostly describes the 



surface brightness profile in the outer region, and R C 2 = 15.5^q '\h 7 ^ kpc and ft = 0.59l 001 



^+0.01 

for another, which describes the central excess brightness. 



5.2. Gas and Gravitating Mass Distributions 

The distribution of the diffuse X-ray emission in different energy bands is roughly sym- 
metric outside to about 4.2' (~ 800/ifg 1 kpc) in MS 0839.9+2938, indicating that the cluster 
is nearly relaxed. In this case the X-ray imaging spectroscopy can be considered as a reliable 
tool to deduce the mass distributions in the cluster. Assuming that the gas is ideal and the 
cluster is in the hydrostatic equilibrium state, and utilizing the best-fit temperature gradi- 
ent obtained in the deprojected analyses and the gas density profile deduced with the two-/3 
model, we calculate the three-dimensional gravitating mass distribution by using 

M(m- kTR ( dlnn 3 i dlnT \ (r\ 
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where M(R) is the total mass within the radius R, G is the gravitational constant, and m p is 
the proton mass. Here the effects of non-thermal turbulent, magnetic pressure, and cosmic 
ray pressure are not taken into account. The resulting distribution of the total gravitating 
mass is shown in Figure 6, along with the 90% errors determined by performing the Monte- 
Carlo simulations that account for the range of temperature and gas profiles allowed by the 
data. The total mass distribution shows an excess beyond that predicted by the (5 model in 
R<30hjQ kpc, which is roughly the optical boundary of the cD galaxy (Nesci et al. 1989) 
and the region where the excess X-ray emission, temperature drop, and iron abundance 
increase are observed. The excess mass (~6x 10 n Mo) is comparable to the mass of the cD 
galaxy as is inferred from its optical luminosity (~lx 1O 12 M within 30/ifo 1 kpc; Nesci et 
al. 1989). At the virial radius -R200 (1640/ifg 1 kpc), within which the average mass density 
is 200 times the current critical density of the universe, the extrapolated gravitating mass 
is M 2 oo = 6.1 x 1O 14 M . In Figure 6, we also show the calculated gas mass fraction along 
with the 90% errors. It increases from about 0.05 at the center to about 0.1 at ~ IOO/i^q 1 
kpc, and then keeps as a constant in the outer regions. For Q^h 2 = 0.024 ±0.001 (Spergel et 
al. 2003), we obtain Q, m = 0.421°;^, which is self-consistent with the value adopted in this 
work. 



With the gas density n g (R) deduced with the best-fit two-(3 model (§5.1), we estimate 
the gas cooling time as 



(cf, Sarazin 1986), where both continuum and line emissions are included (Fig. 6). We find 
that t coo \ ~ 10 8 yr at the center, which is much shorter than the age of the universe at the 
cluster's redshift. The cooling radius i? coo i, defined as the radius where the cooling time is 
equal to the age of the universe, is ~ 150/ifo 1 kpc. This is roughly the radius within which 
the gas temperature begins to descend inwards. 

In order to determine the mass deposit rate of the cooling processes, we fit the spectrum 
extracted in the two innermost annuli (< 7A.lhjQ kpc) where the gas temperature is signif- 
icantly lower than that of the outer regions. We use the cooling flow model MKCFLOW, 
which is absorbed by a column density fixed to the Galactic value (Fig. 7). We find that 
when the low temperature kT\ ow of MKCFLOW is fixed to 0.08 keV, the model yields a 
rather poor fit to the data (xl — 1-69), giving rise to too much cool gas and thus overesti- 
mating the data in <0.8 keV. When /cT low is left free, the fit becomes acceptable (xl = 1.00), 



5.3. Cooling of Gas and Mass Deposit Rate 





-13- 



but it requires kT\ ow = 1. M+o^ keV, which may suggest that the spectrum is dominated 
by a hot component. We then add an additional ambient isothermal component (MEKAL) 
into the model, tying its temperature to kT^^ of MKCFLOW. When /cXi ow is fixed to 0.08 
keV, we obtain kT^^ = 3.40^Q25 keV and M = 53^ M yr~ x ; when kT\ ovl is free, we obtain 
A;T low = 0.56^;| keV, A;T high = 3.42^;^ keV, and M = 67+|M yr" 1 . Whether or not A;T low 
is fixed, the MKCFLOW+MEKAL fit is acceptable (xl = 0.99 and 1.00, respectively), but 
kT low and M are far from being well constrained. On the other hand, the two-temperature 
APEC model also gives a good fit (xl = 0.99) to the data. The derived best-fit low and high 
temperatures are kT low = 2.17^1 keV and kT high = 4.85^1 keV. Thus, given the current 
data quality and instrumental capabilities, the two-phase isothermal model cannot be distin- 
guished from the single MKCFLOW model with a free kT low , or the MEKAL+MKCFLOW 
model. Although not well constrained, the cooling flow model cannot be ruled out for this 
cluster. 



6. Discussions 

6.1. Metal Enrichment in ICM and Role of SNe la in Inner Region 

The inward iron abundance increase found in many nearby clusters and groups around/in 
the giant, dominating elliptical galaxy at the center seems to follow a somewhat universal 
profile. Actually, when we compare the iron abundances in the inner (^70/ifQ 1 ) and outer 
regions (Z Fcin and Z FefiUt , respectively) of the clusters in which an inner abundance en- 
hancement is observed (Buote et al. 2003; Nevalainen et al. 2001; Pratt & Arnaud, 2003; 
Takahashi & Yamashita 2003; Finoguenov et al. 2004; Tamura et al. 2004), we find that 
within the errors the ratio Z-pe^/ Zp ey0nt is roughly a constant as a function of either the ICM 
temperature or the cluster's B-band luminosity (Fig. 8). All these results strongly imply 
that the observed central iron excess has probably been caused by the same enrichment 
process associated with the central dominating galaxy, which is a cD galaxy in most cases. 
Meanwhile, the high spatial resolution observations of the nearby giant ellipticals M87 (Mat- 
sushita et al. 2003) and NGC 5044 (Buote et al. 2003) revealed that the spatial distribution 
of oxygen abundance is approximately uniform out to about 40-70 kpc, inferring that the 
metals synthesized in SNe II are evenly distributed in the cluster /group. Thus, the central 
iron abundance increase should be ascribed mainly to SNe la. 

By performing both the projected and deprojected spectral analyses, we find that in MS 
0839.9+2938 the measured average iron abundance is about 0.4 solar in the outer regions, 
while it increases to ~1 solar within the central 37h^ kpc at the 90% confidence level. 
Reports about detections of such a significant central iron excess in clusters at -2^0.2 are 
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rare in literatures. In 4C+55.16 (z = 0.240; Iwasawa et al. 2001) the metal abundance 
increases from <0.5 solar in outer regions to > 1 solar within the central 34h^Q kpc. In 
Abell 2390 [z = 0.23; Allen et al. 2001), it increases from 0.23^;^ solar in 143 - 7Uh^ 
kpc to 0.48lo'}o solar within < 71h^Q kpc. However, note that the results of the latter two 
cases were acquired by performing projected analysis at the 68% confidence level, and the 
abundance is an average value for both iron and a-elements, although it is iron-dominated. 

In MS 0839.9+2938, the excess iron mass within R < 74hjQ kpc is estimated to be 
M Feiexccss ~ 3.0 x 1O 8 M , which is similar to that of Abell 2390 (M Fe)ex cess ^ 2.5 x 1O 8 M 
within R < 71/ifg 1 kpc, based on the data of Allen et al. 2001). The derived Mpe.excess 
obviously exceeds the iron enrichment caused only by the stellar winds and ram pressure 
stripping of metal-enriched gas (e.g., Trager et al. 2000; Ciotti et al. 1991). Therefore, 
we speculate that a significant part of the excess iron may have been blown into the ICM 
directly by the SN la explosions. 

By reproducing the observed iron-to-silicon abundance ratio in the central 7Ah^ kpc 
of MS 0839.9+2938 (Fe/Si=1.29 + o;^ solar at the 68% confidence level), where a clear iron 
abundance increase is seen, we have attempted to estimate the SN la contribution to the iron 
enrichment. The current supernova models give very different yields to the metals (Gibson 
et al. 1997). When we adopt the theoretical SN la yields obtained with the W7 model and 
the weighted SN II yields in Nomoto et al. (1997), where an IMF slope of x = 1.35 (Salpeter 
1955) and a stellar mass range of 10 — 5OM were assumed for the SN II progenitors, we 
find that the SN la contribution to the total iron mass is 90^}g%, corresponding to an iron 
abundance of ^Fe,SNia = 0.56^o; < j > i solar. This fraction is consistent with those found in NGC 
5044 (67 - 79%; Buote et al. 2003), M87 (55 - 97%; Matsushita et al. 2003), NGC 1399 
(~ 80%; Buote 2002), and the RGH 80 group (~ 80%; Xue et al. 2004). We are not able 
to derive reasonable results by using the SN la yields obtained with either the WDD1 or 
WDD2 model in Nomoto et al. (1997), because both models give rise to too much silicon. 

Using the same method as above we estimate that for the 74 < R < AAAhjQ kpc region, 
where the observed Fe/Si ratio is 0.84^^7 solar (68% confidence), the SN la contributions 
is 75^5%, or ^Fc.sNia = 0.23tg]i7 (W7 model). By comparing these results with those for 
the inner region, we find that the spatial distribution of the SN II contribution to iron is 
roughly constant at about 0.08 solar, which is close to that in Finoguenov et al. (2002; 
^Fe,SNii = 0.1 — 0.15 solar) and Matsushita et al. (2003; Zf c ,snii = 0.2 — 0.4 solar). 

It is crucial to examine if there has been enough time for SNe la to enrich the cluster's 
central region with the observed amount of excess iron. Since the SN la enrichment of iron 
in the ICM is determined both by iron blown out directly into the ICM during the SN la 
explosions and by iron lost in the stellar winds, we calculate the enrichment time t C nrich 
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(Bohringer et al. 2004) using 

^enrich = -^Fe,SNIa/[-^B X (-^SNIa + -Rwind)], (10) 

where Mpe.sNia is the iron mass contributed by SNe la, Lb (2.2 x 10 n /i7 2 L Q ; Nesci et al. 
1989) is cluster's B-band luminosity, -RsNia is the direct SN la iron-enriching rate in units of 
M Q yr^ 1 Lq 1 , and -R W md is the iron-enriching rate due to the stellar mass loss in units of M Q 
yr _1 Lq 1 . As is shown above, we assume that SNe II evenly enriched the gas throughout 
cluster, yielding an iron abundance of 0.08 solar, so by subtracting the SN II contribution, 
we obtain M FCjSNIa = 5.2 ± 0.8 x 1O 8 M . In SN la explosions the rate of iron mass that has 
been directly transferred into the ICM is 

R sma = SR10- 12 L^ &VFe , (11) 

where 77 Fe is the iron yield per SN la event, and SR is the SN la rate in units of SNu (1 
SNu = 1 supernova (1O 1o L_b )0 ) _1 century^ 1 ). In calculations, we choose 7/Fe = 0.5 and 0.7 
M Q , respectively, and adopt either a temporally constant SR of 0.18/iy SNu (Cappellaro et 
al. 1999), or an observation-constrained time-dependent SR that increases with the redshift 
until z ~ 1 and then drops (Dahlen et al. 2004). In the stellar winds, the iron loss rate is 

R wmd = 1.5 x lO^L-y^V, (12) 

where 7F e is the iron mass fraction in stellar winds, and ti$ is the age in units of 15 Gyr 
(Ciotti et al. 1991). We list the calculated enrichment times and redshifts at which the SN 
la enrichment is expected to start in Table 5. We find that when the supernova rate was 
a constant in the past, the derived enrichment time is too close to or even longer than the 
Hubble time at the cluster's redshift (11.1 x 10 9 yr). When the observed time-dependent 
supernova rate provided in Dahlen et al. (2004) is adopted, the derived enrichment times 
are 7.9 Gyr and 6.4 Gyr for rj = 0.5 and 0.7 M Q , respectively. For comparison, we also 
perform similar calculations for the inner 71/ifo 1 kpc region of Abell 2390 and find that the 
resulting enrichment time is 6.6-10.7 Gyr. These results agree with those deduced for the 
Virgo cluster, the Centaurus cluster, the Perseus cluster and Abell 1795 (i C nrich = 5 — 11 
Gyr; Bohringer et al. 2004), indicating that the excess iron can be entirely produced by the 
central brightest galaxy (De Grandi & Molendi 2001; De Grandi et al. 2004). Since only few 
SNe la are observed at z>1.8, we predict that the most distant clusters showing a significant 
central iron excess should be detected at no farther than z ~ 0.5. 



6.2. Scaling Relations 

The derived X-ray temperature, X-ray luminosity and virial mass at i? 2 oo of MS 0839.9+2938 
are kT = 3.65+^ keV, L x = 4.1±°,; 4 x 10 44 erg s" 1 , and M 200 = 6.1 x 10 14 /if 1 M Q , re- 
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spectively. Within the errors, these values are consistent with the observed luminosity-mass 
relation (e.g., Markevitch et al. 1998) and mass-temperature relation (e.g., Sanderson et al. 
2003) after the effect of gas cooling is compensated. Using the gas density obtained with the 
best-fit two-/? model (§5.1) and the analytical power-law temperature profile (§4.4), we cal- 
culate the azimuthally averaged gas entropy S = kT/n 2 ^ 3 in the cluster and plot its spatial 
distribution in Figure 9. In R > 0.03^200 (49/ifo 1 kpc), the profile can be approximated by a 
power-law form with an index of 1.1, as has been expected from shock heating that occurred 
during the spherical collapse (Tozzi & Norman 2001). At 0.1_R 2 oo, the derived gas entropy 
is ~ 167tj|/i7 1/3 keV cm , which is consistent with the entropy-temperature relation for a 
sample of 66 relaxed clusters presented in Ponman et al. (2003). 



7. Summaries 

We present the Chandra ACIS study of the intermediately distant (z = 0.194) cluster of 
galaxies MS 0839.9+2938. By performing both the projected and deprojected spectral anal- 
yses, we find that the emission-measure weighted gas temperature is approximately constant 
at about 4 keV in > 130/iyo 1 kpc. In the inner regions, the gas temperature descends towards 
the center, reaching <3 keV within the innermost 37/ifg 1 kpc. Along with the temperature 
drop, we detect a significant inward iron abundance increase from about 0.4 solar in the 
outer regions to ~1 solar within the central 37hjQ kpc. Besides MS 0839.9+2938, another 
cluster at z ~ 0.2 that has been found to show a similar, strong central iron excess is Abell 
2390 (z = 0.23; Allen et al. 2001). We argue that most of the excess iron is contributed by 
SNe la. By using the observed SN la rate and stellar mass loss rate, we estimate that in 
MS 0839.9+2938 the time needed to enrich the central region with excess iron is 6.4 — 7.9 
Gyr, which is similar to those found for the nearby clusters. In almost the same region 
where the temperature drop and abundance increase are seen, the observed X-ray surface 
brightness profile shows an excess beyond the distribution expected by either the j3 model 
or the NFW model, and can be well fitted with an empirical two-/? model. The origins of all 
these phenomena can be correlated with each other (cf, Makishima et al. 2001). The excess 
amount of gas, a part of which comes from the supernova explosions with heavy metals, 
is trapped around the central dominating galaxy to form the observed surface brightness 
excess. The relatively high gas density and high metallicity both speed up the gas cooling 
process, causing the observed inward gas temperature drop. 

This work was supported by the National Science Foundation of China (Grant No. 
10273009 & 10233040), the Ministry of Science and Technology of China (Grant No. NKBRSF 
G19990754), and Shanghai Key Projects in Basic Research (No. 04JC14079). 
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Table 1. Fits to the Total Spectrum 
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d - 88 -0.26 


44+0-09 
-0.08 

46+° 10 








107.6/91 
102.8/90 


VAPEC 


4.04 (Fixed) 


3 56 +U ' 2i 

q oq + 0.26 
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1 20+ ' 48 
l-42l^ 


42 +,) ' 2y 
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(High) - 
(Low) 
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37+ l) ' 2,/ 
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87.7/87 



Average abundance of the a-elements except for Si. 
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Table 2. Projected and Deprojected Spectral Analyses 



APEC Model 



Region 


Radius 




Projected 






Deprojected 








kT 
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xi 


kT 


Z 
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1 

2 
3 
4 
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kT 
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(keV) 
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0.0-37.3 
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43.9/42 
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1 p=+1.8U 
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34.8/42 
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^' Oy -0.45 


< 0.32 


42+ - 21 
U ' 4Z -0.20 


37.2/37 
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^■' 3 ' -0.57 


< 0.32 
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41.3/42 


5 


230.0-444.1 
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x - '^-l.OO 


qcr+0.65 


30.1/41 











1 Absorptions fixed to the Galactic value 7V H = 4.04 X 10 20 cm 2 (Dickey & Lockman 1990). 



Table 3. Fits to the Spectra Extracted in 0.0-74. lh^ kpc and 74. 1-444. kpc 
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74.1-444.1 


4.04 (Fixed) 
4 7Q+1-28 
4 -' y -1.24 

4.04 (Fixed) 

4 05+ 1 ' 48 
4 ' U0 -1.45 


3-20±g;ig 
3-14±g:il 

q s1 +0.32 
3.81_ 34 

q 01 +0.42 
3.81_ 44 


0.87l°;g 

0.891S-.S 
1 17 +o.73 

1 - ± '-0.60 
1 17+1.73 
J - 1 '-0.60 


U - 4S -0.35 

41 + - 38 
u - 4i -0.34 

+0.44 
u -°'-0.37 
q 7 +0.48 

-0.37 


0.62l°;« 

omt° it 

o.3il8:it 

0-3l±g;ii 


101.2/96 
100.2/95 
79.1/98 
79.1/97 



a Average abundance of the a-elements except for Si. 



Table 4. Fits to the Observed X-ray Surface Brightness Profile 



13 Model 



NFW Model 



i g (10 
37.5 



cm 
+ 1.5 
1.2 



T ) Rc (/ifo 1 kpc) 



l 70 
47 O+0.3 



0.612 



+0.003 
0.003 



101.1/62 



R B (h 70 kpc) 
I95.ll°; 5 5 



6.5 



c 

+0.1 
0.1 



xi 

69.2/62 



Two-/3 Model 

n g ,i (IO" 3 cm" 3 ) Rd (/i^, 1 kpc) ft n g . 2 (IO" 3 cm" 3 ) R c2 (/if,, 1 kpc) & 

21.7l°; 7 65.9l°j 0.660l°;°8s 40.4±i;| 15.5l°; 3 0.5861°;°°? 46.1/59 
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Table 5. Iron Enrichment by SNe la 



SN la Rate 


^enrich(10 9 jr) 


2 a 


(h?o SNu) 


»7 Fo = 0.5 M 


0.7 M 


■q Fc = 0.5 M 0.7 M 


0.18 




q 91 +2-36 


> 1.66 


Observed b 


> 6.78 


6 36+ U yy 


> 1.45 1.30+£« 



a Redshift at which SN la iron enrichment is expected to start. 
b Data are from Dahlcn ct al. (2004). 
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Figure Captions 

Fig.l - (a): Central 4.9' (944/^ kpc) of MS 0839.9+2938 in 0.7-8 keV, which is plotted in 
the logarithmic scale. The image has been background-subtracted, exposure-corrected, and 
smoothed with a minimum significance of 3 and a maximum significance of 5. X-ray intensity 
contours are spaced from 1.4 x 10~ 10 to 4.7 x 10~ 7 photons cm~ 2 s _1 pixel -2 . The inner 
1' region, which is marked by a box, is magnified in (b), where the contour levels are from 
5.7 x 10 -9 to 2.5 x 10~ 7 photons cm~ 2 s _1 pixel -2 , (c): Corresponding DSS optical image 
for the same sky field as in (a) on which the X-ray contours are overlaid for comparison. 

Fig. 2 - Azimuthally-averaged radial surface brightness profile of MS 0839.9+2938 in 0.7-8 
keV. The background is assumed to be spatially uniform and is represented by the dotted 
line. In 15.8 — 800/ifQ 1 kpc, the surface brightness profile can be well described by the /? 
model (dashed line), showing an obvious central emission excess. The fits of the two-/?, /? 
and NFW models to the whole cluster region, and the fit of the /? model to the outer region 
only are shown with the solid, dashed, dot-dashed and dot-dot-dot-dashed lines, respectively. 

Fig. 3 - Two-dimensional distribution of the hardness-ratio of the diffuse emission in MS 
0839.9+2938, which is defined as the ratio of the background-subtracted and exposure- 
corrected counts in 1.4-8.0 keV to those in 0.7-1.4 keV. The cross represents the peak 
of X-ray emission, and the circle (r = 74hjQ kpc) marks the region in which significant 
temperature drop and abundance increase are detected. The image is plotted in linear scale, 
and has been smoothed in the same way as in Figure 1. 

Fig. 4 - (a): Projected (solid) and deprojected (dashed) radial distributions of the gas tem- 
perature, which are obtained with the absorbed APEC model. The best-fit temperature 
profile for the deprojected distribution is shown as a dotted line, (b): The same as (a) 
but for the absorbed VAPEC model, (c): Projected (solid) and deprojected (dashed) radial 
distributions of the average metal abundance, which are obtained with the absorbed APEC 
model, (d): The same as (c) but for the iron abundance obtained with the absorbed VAPEC 
model. 

Fig. 5 - Confidence contours at the 68.3%, 90% and 99% confidence levels for the gas tem- 
perature and iron abundance derived with the deprojected VAPEC model for the innermost 
(R < 37.3/ifo kpc) and the third (74.1 < R < 130.1/ifo kpc) annular regions. 

Fig. 6 - Upper and middle panels: calculated radial distributions of the total gravitating 
mass, gas mass and gas mass fraction, along with the 90% error limits (shown as dotted 
lines). The total mass was calculated by using both /? and the best-fit two-/? model. The gas 
mass and gas fraction were calculated by using the best-fit two-/? model. Lower panel: the 
calculated cooling time as a function of radius. The horizontal, dashed line represents the 
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age of the universe at the redshift of cluster. Radii within which the gas temperature drops 
and iron abundance increases significantly (i? a ), and at which the gas temperature begins to 
decrease inwards (i?b), are indicated by vertical lines. 

Fig. 7 - Upper and middle panels: Spectrum extracted in < 74.1/ifQ 1 kpc, where the gas 
temperature is significantly lower than that of the outer regions, along with the best-fit two- 
temperature APEC model and residuals. Lower panel: model ratios of the single MKCFLOW 
model and MEKAL+MKCFLOW model to the two-temperature APEC model. kT low of the 
MKCFLOW component is fixed to 0.08 keV. 

Fig. 8 - Ratios of the iron abundance in cluster's inner region to that in the outer region 
Zin/Z ou t as a function of the cluster's (a) average gas temperature and (b) optical luminosity 
for nearby clusters (solid; Tamura et al 2004; Girardi et al. 2000; David et al. 1995; Pratt & 
Arnaud 2003) and two intermediately distant clusters, i.e., MS 0839.9+2938 (dashed) and 
Abell 2390 (dot-dashed. This is the average abundance ratio rather than iron abundance 
ratio; Allen et al. 2001). The ratio for the nearby clusters is shown as a horizontal dotted 
line. 

Fig. 9 - Radial distribution of the gas entropy (crosses and the dashed line) along with the 
theoretical expectation quoted from Tozzi et al. 2001 (solid line). 
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Fig. 1.— (a): Central 4.9' (944/^ kpc) of MS 0839.9+2938 in 0.7-8 keV, which is plotted in 
the logarithmic scale. The image has been background-subtracted, exposure-corrected, and 
smoothed with a minimum significance of 3 and a maximum significance of 5. X-ray intensity 
contours are spaced from 1.4 x 10~ 10 to 4.7 x 10~ 7 photons cm -2 s -1 pixel -2 . The inner 
1' region, which is marked by a box, is magnified in (b), where the contour levels are from 
5.7 x 10~ 9 to 2.5 x 10~ 7 photons cur 2 s -1 pixel -2 , (c): Corresponding DSS optical image 
for the same sky field as in (a) on which the X-ray contours are overlaid for comparison. 
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Fig. 2. — Azimuthally-averaged radial surface brightness profile of MS 0839.9+2938 in 0.7-8 
keV. The background is assumed to be spatially uniform and is represented by the dotted 
line. In 15.8 — 800/ifQ 1 kpc, the surface brightness profile can be well described by the (3 
model (dashed line), showing an obvious central emission excess. The fits of the two-/?, (3 
and NFW models to the whole cluster region, and the fit of the f3 model to the outer region 
only are shown with the solid, dashed, dot-dashed and dot-dot-dot-dashed lines, respectively. 
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Fig. 3. — Two-dimensional distribution of the hardness-ratio of the diffuse emission in MS 
0839.9+2938, which is defined as the ratio of the background-subtracted and exposure- 
corrected counts in 1.4-8.0 keV to those in 0.7-1.4 keV. The cross represents the peak 
of X-ray emission, and the circle (r = 74hjQ kpc) marks the region in which significant 
temperature drop and abundance increase are detected. The image is plotted in linear scale, 
and has been smoothed in the same way as in Figure 1. 
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Fig. 4. — (a): Projected (solid) and deprojected (dashed) radial distributions of the gas 
temperature, which are obtained with the absorbed APEC model. The best-fit temperature 
profile for the deprojected distribution is shown as a dotted line, (b): The same as (a) 
but for the absorbed VAPEC model, (c): Projected (solid) and deprojected (dashed) radial 
distributions of the average metal abundance, which are obtained with the absorbed APEC 
model, (d): The same as (c) but for the iron abundance obtained with the absorbed VAPEC 
model. 




Gas Temperature (keV) 

Fig. 5. — Confidence contours at the 68.3%, 90% and 99% confidence levels for the gas tem- 
perature and iron abundance derived with the deprojected VAPEC model for the innermost 
(R < 37.3/ifQ 1 kpc) and the third (74.1 < R < 130. l/i^ 1 kpc) annular regions. 



-27- 




Fig. 6. — Upper and middle panels: calculated radial distributions of the total gravitating 
mass, gas mass and gas mass fraction, along with the 90% error limits (shown as dotted 
lines). The total mass was calculated by using both (3 and the best-fit two-/? model. The gas 
mass and gas fraction were calculated by using the best-fit two-/? model. Lower panel: the 
calculated cooling time as a function of radius. The horizontal, dashed line represents the 
age of the universe at the redshift of cluster. Radii within which the gas temperature drops 
and iron abundance increases significantly (R a ), and at which the gas temperature begins to 
decrease inwards (Rb), are indicated by vertical lines. 




Fig. 7. — Upper and middle panels: Spectrum extracted in < 74.1/i^q 1 kpc, where the gas 
temperature is significantly lower than that of the outer regions, along with the best-fit two- 
temperature APEC model and residuals. Lower panel: model ratios of the single MKCFLOW 
model and MEKAL+MKCFLOW model to the two-temperature APEC model. kT low of the 
MKCFLOW component is fixed to 0.08 keV. 
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Fig. 8. — Ratios of the iron abundance in cluster's inner region to that in the outer region 
Zm/Zout as a function of the cluster's (a) average gas temperature and (b) optical luminosity 
for nearby clusters (solid; Tamura et al 2004; Girardi et al. 2000; David et al. 1995; Pratt & 
Arnaud 2003) and two intermediately distant clusters, i.e., MS 0839.9+2938 (dashed) and 
Abell 2390 (dot-dashed. This is the average abundance ratio rather than iron abundance 
ratio; Allen et al. 2001). The mean ratio for the nearby clusters is shown as a horizontal 
dotted line. 
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Fig. 9. — Radial distribution of the gas entropy (crosses and the dashed line) along with the 
theoretical expectation quoted from Tozzi et al. 2001 (solid line). 



